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1 Introduction 



The angle (3 of the unitarity triangle has been determined to sin(2/3) = 0.725 ± 0.037 
[1] from time-dependent CP asymmetries in b — > ccs transitions. If sub-leading de- 
cay amplitudes can be neglected as argued in [2], time- dependent CP asymmetries in 
penguin-dominated b — > sqq transitions should also take the value ± sin(2/3). There now 
exist various measurements [3], which on average point to the significantly smaller value 
0.43 ±0.07. It is not inconceivable that flavour-specific new flavour- violating interactions 
cause anomalous effects in b — > s transitions without resulting in inconsistencies with 
other measurements. This would be a rather spectacular resolution of the apparent dis- 
crepancy. But before this conclusion can be drawn, a thorough study of the sub-leading 
decay amplitudes is necessary to ascertain the Standard Model expectation. This is 
undertaken here in the framework of QCD factorization [4]. 

The analysis is based on the next-to-leading order (NLO) factorization calculations 
performed in [5] , where numerical values of the time-dependent CP asymmetries for the 
4>K$ and rj'Kg final states have already been given. In this Letter I include a larger 
set of final states (see also the recent work [6, 7]), and consider a more detailed error 
estimation that includes a scan of the theoretical parameter space [8]. I also discuss 
constraints on the sub-leading decay amplitudes that do not rely on factorization but 
are inspired by it. Another method to constrain the differences of time-dependent CP 
asymmetries in b — > ccs and b — > ssq transitions based on systematic approximations to 
the strong interactions relies on the assumption of SU(3) flavour symmetry. This results 
in bounds on the magnitude of this difference, but the sign cannot be determined [9]. 
An estimate in a model of (long-distance?) final state interactions is given in [6]. 

The time-dependent CP asymmetry in decays to CP eigenstates is given by 



Br(B°(t)^ /) - Br(i? (t) ^ f) 
Br(S°(t)-/) + Br(fl°(t)-/) 



Sf sin(Ams t) — Cf cos(Ajtib t), 



(1) 



with Attib the B B mass difference. The B decay amplitude involves two weak cou- 
plings VpbV* and two strong interaction amplitudes a p f . I write 



A(B - /) = V cb V c * s a) + V ub V: s a} cx 1 + e"*? d 



where 



df = 6km— = tKNidf with ckm 
a f 



VcbV* 



A standard calculation now gives 



0.025. 



(2) 
(3) 



AS f = -rjfSf - sin(2^) 



2 Re(d f ) cos(2/3) sin 7 + \d f \ 2 (sin(2/3 + 2 7 ) - sin(2/?)) 



.4 



CP,/ 



1 + 2 Re(df) cos7 + \df\ 2 

2 lm(df) sin 7 

2 Ke(df) cos7 + |<i/| 2 



(4) 



(5) 



f 



Here rjf denotes the CP eigenvalue of /. (All final states discussed below have rjf = — 1.) 
The quantity ASf is the central object of this Letter. One notes that (a) df is suppressed 
by a small ratio of CKM elements, €km, leading to the expectation that —r]fSf ~ sin (2/3) 
(see above); (b) if df is small as expected, then to first order in df the two asymmetries Sf 
and Cf involve independent hadronic parameters, namely the dispersive and absorptive 
part of df = a U fja C f. 



2 Anatomy of AS/ in factorization 

The hadronic amplitudes a p f, p = u,c are sums of "topological" amplitudes, referring to 
tree (T, C), QCD penguin (P p ), singlet penguin (S p ), electroweak penguin (P^w ^ewc) 
and annihilation contributions. The relation to the "flavour" amplitudes used in QCD 
factorization [5] is T <-> a±, C <-> a 2 , -P p <-> a^ + ^f, <-> a3 + /3g 3 , and (-Pgw ^ew c*) ^ 
(^3 EW a 4Ew) with the difference that the exclude form factors, decay constants and 
the CKM factors, while the topological amplitudes exclude only the CKM factor. In 
addition, a penguin amplitude such as P c may be a sum of several a p terms depending 
on the flavour flow to the final state. The expressions for all relevant decay amplitudes 
in terms of flavour amplitudes are collected in Appendix A of [5] . Schematically, for the 
strangeness-changing decays B° —>■ MK°, the hadronic amplitude ratio is given by 

df ~ £km pc , (oj 

where the dominant amplitudes have been indicated. Note that the amplitudes P u , C, . . . 
depend on the final state /. 

In the QCD factorization framework the topological amplitudes are computed in the 
form [4] 

T, C, P c > u , ... = £ C(ji h ) x {F BM ^ x T\fi h , fi 8 ) *f M ,$ M M 

terms " 

l+a s +... 

+ Jb®B (// s ) * T U {fih,^l\ * J n {fil,fJLs) * fhh ®Mi (Ps) * /M 2 $M 2 (^*) } 

a s +... 

+ 1 /m b -suppressed terms (7) 

reducing the hadronic input to form factors F BM and light-cone distribution amplitudes 
$x- The underbraces indicate the order in perturbation theory to which the various 
short- distance kernels are computed at NLO. The numerical implementation of (7) also 
includes some l/m& power corrections from scalar penguin operators, and from an esti- 
mate of annihilation topologies. The accuracy of the treatment is generically limited by 
Aqcd/^?) ~ (10 — 20)% at the amplitude level. 

The actual uncertainties affect different observables to a different degree and must be 
estimated on a case-by-case basis. The "colour-allowed" amplitudes T, -P^w are rather 
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certain, while the "colour-suppressed" amplitudes C, Pew c receive contributions from 
spectator scattering (the second line of (7)) enhanced by large Wilson coefficients, and 
are inflicted by larger uncertainties. The QCD penguin amplitudes include uncertain 
annihilation contributions, although the ratio P u / P c is less affected. Finally, the singlet 
amplitude S p involves several specific decay mechanisms [10], which are difficult to com- 
pute quantitatively, though none of them seems to be of particular importance for the 
CP asymmetries. Eq. (6) indicates that ASf involves some of the less certain amplitudes. 

The numerical analysis below takes into account all flavour amplitudes following [5], 
but it suffices to focus on a few dominant terms to understand the qualitative features 
of the result. Then, for the various final states, the relevant hadronic amplitude ratio is 
given by 

71 Ks d f — p Ks d f rp*j — 

V'K S d f ~ tn__M ^ Ks (8) 

V A H + ^ K A H + W 

T]K S d f ujK s df 

The convention here is that quantities in square brackets have positive real part. (Recall 
from (4) that ASf mainly requires the real part of df.) In factorization Ke[P u /P c ] 
is near unity, roughly independent of the particular final state, hence ASf receives a 
nearly universal, small and positive contribution of about 2€km cos(2/5) sin 7 « 0.03. On 
the contrary the magnitudes and signs of the penguin amplitudes' real parts can be very 
different. Ignoring uncertainties, I find |Re[P c ]| in the proportions 

7r°K : p°K : rj'K : <j>K : rjK : uoK 
1 : 0.5 : 2.2 : 0.8 : 0.5 : 0.5 

Hence the influence of the colour-suppressed tree amplitude C determines the difference 
in ASf between the different modes. For (jr°, 77, oS)Ks the effect of C is constructive, but 
for (p, rj')Ks it is destructive. However, the magnitude of Re [P c ] is much larger for rj'Kg 
than for pK$, hence lie (df) remains small and positive for the former final state, but 
becomes negative for the latter. 



3 Factorization results 

The result of the calculation of ASf is shown in Table 1. The column labeled u ASf 
(Theory)" uses the input parameters (CKM parameters, strong coupling, quark masses, 
form factors, decay constants, moments of light-cone distribution amplitudes) summa- 
rized in Table 1 of [5]. In particular \V ub /V cb \ = 0.09 ± 0.02 and 7 = (70 ± 20)° is used. 
The uncertainty estimate is computed by adding in quadrature the individual parameter 
uncertainties. The central values are in good agreement with those given in [6], which 
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Mode 


ASf 


(Theory) 


ASf [Range] 


Experiment [3] (BaBar/Belle) 


n°K s 





nv+o.05 
yj 1 —0.04 


[+0.02,0.15] 


-0.39i°;f? 


(-0.38±8:g/-0.43±gS) 


P°K S 


-0 


00+O.O8 
U °-0.12 


[-0.29,0.02] 






v'Ks 





r, -,+0.01 


[+0.00,0.03] 


-o.3o±8:g 


(-0.43±8;^/-0.07±g;iI) 


T]Kg 





in +o.n 

iu -0.07 


[-1.67,0.27] 






<PK S 





no+o.oi 

uz -0.01 


[+0.01,0.05] 




C-o.23±8ii/-o.67±8:iS) 


uK s 





i q+0.08 
io -0.08 


[+0.01,0.21] 


-0.18i°;| 


(-0.23i°;^/+0.02l^) 



Table 1: Comparison of theoretical and experimental results for ASf. 

also uses the input from [5]. For the final states p°Kg and uKg they differ from those 
given in [7], where the leading order (naive factorization) approximation is employed, 
and the electroweak penguin amplitudes are neglected. The next-to-leading order cor- 
rection included in the present calculation has a large impact on the branching fractions 
of penguin-dominated modes and is crucial for a successful comparison of QCD factor- 
ization results with data. Nonetheless, the NLO correction to ASf is never larger than 
about 30%, since the amplitude enhancement partially cancels in the ratio df. The NLO 
correction also eliminates the large renormalization scale uncertainty present at leading 
order. 

The result displays the anticipated pattern. The variation of the central value from 
the nearly universal contribution of approximately €km is due to Re [C/P c ], and the error 
comes primarily from this quantity. It is therefore dominated by the uncertainty in the 
hard-spectator scattering contribution to C, and the penguin annihilation contribution to 
P c . In general one expects the prediction of the asymmetry Sf in factorization to be more 
accurate than the prediction of the direct CP asymmetry Cf, since Sf is determined by 
Re (df/a°f) which is large and calculated at next-to-leading order, while Cf is determined 
by Im (af/a,f), which is small and currently known only at leading order. The resultant 
error on ASf is roughly of the size of ASf itself. Since this is small, one arrives at 
accurate constraints, in particular for the final states rj'Kg and 4>Kg. It is striking that 
the theoretical prediction of ASf is positive, with the exception of p°Kg, while the 
experimental data are all negative. 

Quadratic addition of theoretical errors may not always lead to a conservative error 
estimate. Furthermore, the default parameters adopted in [5] do not lead to the best 
description of the data. As shown there, a different choice of a few parameters (defining 
certain "scenarios") results in a very good description of data - however, some observ- 
ables, in particular the colour-suppressed tree amplitude C important to the present 
discussion, then take values outside the range estimated by quadratic error estimation. 
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To allow for this possibility I perform a random scan of the allowed theory parame- 
ter space. For any observable I take the minimal and maximal value attained in this 
scan to define the predicted range of this observable. However, in doing so I discard all 
theoretical parameter sets which give CP-averaged branching fractions not compatible 
within 3 sigma with the experimental data, that is I require 8.5 < 10 6 Br (7r°Ji ) < 14.5, 
0.3 < 10 6 Br (p°K°) < 9.9, 5.3 < 10 6 Br (<f>K°) < 11.9, 2.9 < 10 6 Br (luK°) < 8.3, 
10 6 Br (r]K ) < 6.0. No further condition is imposed, neither from the corresponding 
charged decay modes, nor any other decay, or from direct CP asymmetries (since these 
depend on other hadronic parameters as mentioned above). Note that I also do not 
require the theoretical parameters to reproduce the rj'K branching fraction. The reason 
for this is that in [5] the singlet contribution F 2 to the B — > rj' form factor is set to 
zero simply for lack of better information. Since a non-zero F 2 can affect the branching 
fraction significantly [10], requiring the r]'K branching fraction to reproduce the data 
for F 2 = would be overly restrictive on the remaining theory parameter space. Never- 
theless, one finds that the distribution of B° — > rj'K branching fractions generated by 
the models that survive the other branching fraction restrictions has a (broad) maximum 
at 67 ■ 10~ 6 in nice agreement with experimental data. 

The resulting ranges for ASf from a scan of 200000 theoretical parameter sets is 
shown in the column labeled "AS"/ [Range]" in Table 1. It is seen that the ranges are 
in fact not much different from those obtained by adding parameter uncertainties in 
quadrature - except for the r)Kg final state, for which almost any value of Sf is possible. 
To understand this exception, one must know that similarly large ranges can appear 
also for other final states when no branching fraction restriction is imposed. These large 
values of ASf originate from small regions of the parameter space, where by cancellations 
the leading penguin amplitude P c becomes very small. This leads to large amplifications 
of C/P c , and hence ASf. Such small values of P c always lead to very small branching 
fractions, hence they are excluded by observations except for the case of r/Ks, where no 
lower limit on the branching fraction exists at present. 

The parameter scan contains more interesting pieces of information than the ranges 
of ASf, since it allows to establish correlations between ASf and input parameters, 
between the ASf for different final states etc. in the framework of QCD factorization. For 
instance, one finds that the "good" models prefer a strange quark mass around 80 MeV, 
smaller renormalization scales and a moderate annihilation contribution pA ~ 0.7e lr ^ A 
with \<Pa\ < 70°, all of which affects the magnitude of the dominant QCD penguin 
amplitude. Space does not permit a detailed discussion here, but Figure 1 shows the 
correlation between ASf and the direct CP asymmetry Cf (see (1), (5)) taking / = <f)K° 
and ujK° as examples. The distribution of points (each corresponding to one theoretical 
parameter set) does not reveal any particular correlation between the two observables, 
especially after the branching fraction restriction, as could have been guessed from the 
fact that they mainly involve independent hadronic parameters. The Figure also shows 
that the requirement that the experimental branching fractions be reproduced within 3 
sigma narrows the distribution considerably. Similar conclusions apply to all other final 
states. 
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Figure 1: Correlation between ASf and C/ (direct CP asymmetry) for / = <pK° (left) 
and / = uK° (right). Theory parameter models compatible with the experimental 
branching fractions (as described in the text) are in grey (red), all others in black. 
Based on a sample of 50000 input parameter models. 



4 Discussion 



Given the important role of ASf in the detection of anomalous b — > s flavour transitions, 
one may question the assumptions that go into the factorization approach or attempt 
to find independent validations. Also, given the current experimental status, it would 
already be interesting to know that ASf should be positive, no matter its precise value. 
Can one establish ASf > (except for pK$) with little assumptions on hadronic physics? 
Recall from (4) that ASf is roughly 

2e K MCOs(2/?)sin 7 Re . (10) 

Large enhancements relative to the factorization predictions require an enhancement of 
the hadronic amplitude ratio. The first option is a strong suppression of a C f, but this is 
excluded by the branching fraction measurements (see also the discussion in the previous 
section) . The second option is an enhancement of a J by a factor of several. Can this be 
excluded, or can at least the sign of Re (a^/a'j) be determined? 

The only approach to non-leptonic decays other than factorization based on a small 
expansion parameter uses SU(3) flavour symmetry to relate amplitudes of final states 
belonging to the same SU(3) multiplet. In applications of the method to ASf one uses the 
branching fractions of b — > d transitions to bound \df\ of the related b — > s transitions [9]. 
The best possible limit in this method is \df \ < X 2 ~ 0.05 (A the Wolfenstein parameter), 
so the theoretical limit of this method is \ASf \ <^ 0.07. In practice, depending on the 
values of the b — > d branching fractions and the final state /, the bound is considerably 
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weaker, although the region of interesting values (indicated by the factorization results) 
may eventually be approached for some final states. Note that the sign of ASf is not 
determined by this method. 1 

A limited amount of information can be obtained from final states related to the 
given one by isospin symmetry, or from other observables related to the given final state. 
As already mentioned above, the measurement of the direct CP asymmetry (Cf) is of 
limited use if it is small, since it constrains the imaginary part of cif/a'f rather than 
the real part. On the other hand, a very large direct CP asymmetry (for <f)Ks, rj'Ks, 
ti°Ks) would suggest that Re(d/) could also be large, but this is not rigorous. It would 
certainly imply large violations of factorization, and hence cast doubt on the results in 
Table 1. No such large direct CP asymmetries have been observed to date for the final 
states discussed here. 

The asymmetry Sf is more closely related to ratios of CP-averaged branching frac- 
tions, which also depend mainly on real parts of amplitude ratios. In the following I 
consider the pairs (MK° , MK~), including the charged partners of M for M = it, p. 
The decay amplitudes can be parameterized as 

A{M-R°) = P + e" i7 P u 
V2A(M°K-) = [P + P EW ] + e~ l ~< [T + C + P u ] 

A(M + K~) = [P + P C ' EW ] + e^ 7 [T + P u ] 
V2A(M°K°) = [-P + P EW - P C > EW ] + e~ 17 [C - P u ] (11) 

for M = 7T,p (assuming isospin symmetry), and 

A(MK~) = [P + P C > EW ] + e" 17 [T + C + P u \ 

A(MK°) = P + e~ i7 [C + P u ] (12) 

for the remaining M = rf'\ 0, to. The notation is chosen so that it indicates the dominant 
contribution to each amplitude; the dependence of P, P u , ... on M is not spelled out. 
In (12) the CKM-suppressed penguin amplitude P u is redundant and could be absorbed 

1 It may be noted that the application of the SU(3) approach to final states containing 77, 77', u> and (ft 
requires additional assumptions beyond SU(3). In the SU(3) limit these mesons would be pure octet or 
singlet states, but reality is far from this limit, in particular in the case of w and </>, which are believed 
to be pure up-down and strange quark states, respectively. In [9] this SU(3) breaking singlet-octet 
mixing effect is taken into account by assuming that the operator matrix elements with the physical 
meson states are related to those with the putative SU(3) states by a single mixing angle. This is an 
assumption that cannot be justified in any controlled approximation [11]. Rather one must introduce a 
separate mixing angle for every operator. The existence of large mixing for uj — <f> and 77 — 77' should be 
taken as an indication that a SU(3) treatment might be unreliable, since for every operator a separate, 
presumably large, mixing angle must be introduced. Phenomenological evidence related to the matrix 
elements of current operators may indicate that this SU(3) breaking effect is nearly universal and could 
be described by a single mixing-angle in the quark-flavour-basis [12], but little is known about the matrix 
elements of the effective weak Hamiltonian. 
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Modes 


t 


c 


p" 


p EW 


pC,EW 


nK 


-0.13 


-0.06 


0.02 


0.13 


0.03 


pK 


0.27 


0.13 


0.01 


-0.29 


-0.07 


rj'K 


-0.03 


-0.01 


0.02 




0.01 


rjK 


0.34 


0.14 


0.02 




-0.05 


<j>K 


0.01 


0.00 


0.02 




0.01 


ujK 


0.23 


0.11 


0.02 




-0.08 



Table 2: Estimates of the real part of the amplitude ratios x in scenario IV of [5]. 



into C. For M = <fi the "tree" amplitudes T, C are actually annihilation amplitudes 
and thus very small, provided is a pure ss state, as will be assumed here. It is clear 
from (12) that nothing can be learned from the charged decay for M = 77^, 0, lu without 
additional assumptions, since it involves two new amplitudes (the colour-suppressed elec- 
troweak penguin p c ' EW ; and T). However, I shall now expand the ratios of CP-averaged 
branching fractions under the premise that certain amplitude ratios are small. To this 
end, note that T,C,P U which multiply e _n are proportional to €km, while the elec- 
troweak penguin amplitudes are suppressed by the electromagnetic coupling. Defining 
x = X/P and counting €km ~ A 2 with A a counting parameter of order 1/5, the natural 
magnitudes of the amplitude ratios are t,p EW ~ A, and c, 

p u^C,EW „ 

A 2 . Estimates of 

the real parts of the amplitude ratios are given in Table 2 using the scenario IV of [5] as 
input. In the following discussion, c and p u are allowed to be enhanced to order A. 

Turning first to / = rf^K, 0K,ujK, Eq. (12) implies that even when an enhancement 
of the amplitudes c, p u by a factor of several to order A is allowed, they do not appear in 
the ratio of CP-averaged branching fractions at first order in A. Thus, with an accuracy 
of a few percent, 

t b o Br (M°K~) 



R(f) 



1 + 2cos7Re(t). 



(13) 



r B+ Br (M°K°) 

Hence Re (t) can be determined from data, if R(f) is sufficiently different from 1 (to 
justify the neglect of the order A 2 terms), but ASf oc Re (c+p u ). The colour- allowed tree 
amplitude T is believed to be well-predicted in factorization, and has a small absorptive 
part. Assuming this, an accurate measurement of R(f) for / = tj^'K^ujK provides an 
estimate of Re (P), of which the sign should be reliable. Making the same assumption 
for C constrains the contribution from Re (c) to ASf, but in this case the assumption 
is already questionable. The contribution from Re (p u ) is not constrained as long as 
it is of order A. However, one may argue that if Re (p u ) is enhanced to order A by 
whatever mechanism, then - probably - the absorptive part Im (p u ), and hence the direct 
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CP asymmetry, will also be of order A. Similar arguments can be applied to the ttK and 
pK system (11). To linear order in A, 



The electroweak penguin amplitudes are now important. For pK the corresponding pref- 
actor reduces the branching fraction ratio by a factor of three. In fact, the contribution 
is so large that the linear approximation becomes inapplicable to the pK final state. 
For ttK, the complete set of three branching fraction ratios can be used in principle to 
determine the real parts of t, c and p EW simultaneously with a relative uncertainty of 
order A in the linear approximation. However, the current experimental nK data does 
not lead to useful results. 

I conclude from this discussion that it is very difficult to constrain ASf independent 
of theoretical assumptions using only experimental data (other than the measurement of 
ASf itself). With some plausible dynamical assumptions bounds can be derived using 
SU(3), or the real parts and signs of amplitudes related to the quantities of interest can be 
determined and compared to the factorization calculations, thus providing cross-checks. 

5 Conclusion 

QCD factorization calculations of the time-dependent CP asymmetry in hadronic b — > s 
transitions yield only small corrections to the expectation —rjfSf ~ sin(2/?). With the 
exception of the p°Kg final state the correction ASf is positive, slightly strengthening 
the discrepancy with the current average experimental data. The effect and theoretical 
uncertainty is particularly small for the two final states (j)Ks and rj'Ks already analyzed 
in [5]; the calculation of ASf for the final states p°K$ and rjKs, however, is more 
susceptible to errors because of amplitude cancellations. The final-state dependence of 
ASf is ascribed to the colour-suppressed tree amplitude. 

It appears difficult to constrain ASf theory-independently by other observables. In 
particular, the direct CP asymmetries or the charged decays corresponding to / = MK$ 
probe hadronic quantities other than those relevant to ASf, if these observables take val- 
ues in the expected range. Large deviations from expectations such as large direct CP 
asymmetries would clearly indicate a defect in our understanding of hadronic physics, 
but even then the quantitative implications for Sf would be unclear. A hadronic in- 
terpretation of large ASf would probably involve an unknown long-distance effect that 
discriminates strongly between the up- and charm-penguin amplitude resulting in an en- 
hancement of the up-penguin amplitude. No model is known to me that could plausibly 
produce such an effect. 
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